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Abstract. The preceding paper [J. D. Biggs and J. A. Cina, xxx; referred to here as 
Paper 1], describes a strategy for externally influencing the course of short- time 
electronic excitation transfer (EET) in molecular dimers and observing the process by 
nonlinear wave-packet interferometry (nl-WPI). External influence can, for example, be 
exerted by inducing coherent intramolecular vibration in one of the chromophores prior 
to short-pulse electronic excitation of the other. Within a sample of isotropically oriented 
dimers having a specified internal geometry, a vibrational mode internal to the acceptor 
chromophore can be preferentially driven by electronically nonresonant impulsive 
stimulated Raman (or resonant infrared) excitation with a short polarized "control" pulse. 
A subsequent electronically resonant polarized pump then preferentially excites the 
donor, and EET ensues. Paper 1 investigates control-pulse-influenced nl-WPI as a tool 
for the spectroscopic evaluation of the effect of coherent molecular vibration on 
excitation transfer, presenting general expressions for the nl-WPI difference signal from a 
dimer following the action of a control pulse of arbitrary polarization and shape. 
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Electronic excitation is to be effected and its inter-chromophore transfer monitored by 
resonant pump and probe "pulses," respectively, each consisting of an optical-phase- 
controlled ultrashort pulse-pair having arbitrary polarization, duration, center frequency, 
and other characteristics. Here we test both the control strategy and its spectroscopic 
investigation — with some sacrifice of amplitude-level detail — ^by calculating the pump- 
probe difference signal. That signal is the limiting case of the control-influenced nl-WPI 
signal in which the two pulses in the pump pulse-pair coincide, as do the two pulses in 
the probe pulse-pair. We present calculated pump-probe difference signals for (1) a 
model excitation-transfer complex in which two equal-energy monomers each support 
one moderately Franck-Condon active intramolecular vibration; (2) a simplified model of 
the covalent dimer dithia-anthracenophane, representing its EET dynamics following 
selective impulsive excitation of the weakly Franck-Condon active anthracene 

vibration at 385 cm"^; and (3) a model complex featuring moderate electronic-vibrational 
coupling in which the site energy of the acceptor chromophore is lower than that of the 
donor. 

1. BACKGROUND & SET-UP 

Recent theoretical studies of electronic excitation transfer (EET) have explored 
the influence of coupling between the purely electronic degrees of freedom and 
"environmental" nuclear motion (both intra- and intermolecular). Making use of 
experimentally guided choices for bath-mode correlation functions, Aspuru-Guzik and 
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co-workers ' performed computations tracking useful measures for the influence of 
different environment-induced processes on energy transfer. Plenio et al?''^ subsequently 
analyzed the fundamental effects of electronic dephasing, relaxation, and trapping within 
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the context of models reljdng on a phenomenological Markovian description of the 
environmental interactions. Both groups have made important headway in elucidating the 
observed interplay among coherent electronic dynamics, medium-induced decoherence 
and relaxation, and trapping processes in complex multi-chromophore systems such as 
the Fenna-Matthews-Olson complex, ' bacterial photosynthetic reaction centers, and 
conjugated polymers.^ In the context of these studies, the electronic degrees of freedom 
are regarded as the "system," while the "bath" comprises both inter- and intramolecular 
vibrational degrees of freedom. 

Other studies have addressed certain manifestly non-Markovian effects that can 
occur in the presence of fairly strong coupling between electronic excitations and 
intramolecular or local vibrations. Jang, Jung, and Silbey^ formulated a time-dependent 
generalization of Forster theory accounting for the effects on energy treinsfer of coherent 
vibrational motion in a Franck-Condon active mode. Cina and Fleming^*' performed 
model-system cdculations seeking to rationalize the observation of vibrational coherence 
transfer in time-resolved polarized fluorescence up-conversion measurements on LH-1.^^' 
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In the preceding paper^^ (referred to as Paper 1), Biggs and Cina describe a 
possible means of exerting external influence over the time-course of initial electronic 
excitation transfer by inducing a coherent intramolecular vibration in the acceptor 
chromophore prior to short-pulse excitation of the donor. Their strategy builds on the 
work of Jang, Jung, and Silbey*^ and Cina and Fleming, and seeks to control the 
instantaneous "rate" of excitation transfer by initiating vibrational motion in a Franck- 
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Condon active acceptor mode in order to render time-dependent the acceptor absorption 
frequency relative to the emission frequency of the donor. 

Motivated by current ultrafast experimental capabilities,^'^'^'^'^^'^^'^'^'^^'^^'^"'^^ this 
recent theoretical work generalizes Forster's foundational insights on excitation 
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transfer, and builds on other detailed studies of the dynamics of few- or multi-level 
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quantum mechanical systems in condensed media. • > > • > 

Paper 1 presents general expressions for the nl-WPI difference signal from an 
EET dimer following the action of a control pulse of arbitrary polarization and shape. 
The signal expressions given there apply to a sample of isotropically distributed dimers, 
having a specified internal geometry, in which the acceptor vibration can be 
preferentidly driven by impulsive stimulated Raman excitation with a short, polarized 
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"control" pulse. ' ' An electronically resonant, polarized "pump" comprising two 
optical-phase-related pulses then preferentially excites the donor, and the action of a 
polarized phase-controlled "probe" pulse-pair enables the detection of a wave-packet 
interferometry signal. Such control-pulse-influenced nl-WPI measurements will provide 
highly detdled nuclear-amplitude-sensitive records of the effect of the externally induced 
nuclear motion on the time-course of EET. Here, we perform numerical tests of this 
approach to the control and observation of short-time EET by applying the previous 
general treatment to some simplified model systems, and by collapsing the pump and 
probe pulse-pairs in the WPI sequence each to a single pump or probe pulse. 
Paper 1 treats a molecular dimer with a field-free Hamiltonian 

H=|0>//o(0| + |l)//i(l| + |l'>^r{l1 + |2>//2(2| + /{|l')(l| + |l)(l1} . (1) 
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H^(n = 0,1, 1', 2 ) is the nuclear Hamiltonian in the uncoupled electronic state | n) , 
where 



0) = \gg) 


(both monomers unexcited) 


(2) 




("donor" excited) 


(3) 




("acceptor" excited) 


(4) 


2) = \ee} 


(both monomers excited) , 


(5) 



and / is a parameter determining the strength of the energy-transfer coupling (taken here 
to be constant). The designation of one molecule as the donor and the other as the 
acceptor is arbitrary in the case of a homo-dimer of monomers with equal site energy. 
The signal calculations shown below for such systems treat both chromophores 
identically; their excitation and de-excitation probabilities depend upon the direction 
cosines between their individual transition dipole moments and the various field 
polarizations, hi calculations of time-resolved signals from hetero-dimers, the donor 
(acceptor) will be the monomer having the higher (lower) site energy. 

The Hamiltonian (1) is general in form, allowing any number of intramolecular 
vibrational and bath degrees of freedom — including "mode-sharing" between the 
monomers — ^but we specialize here to some simple cases allowing illustrative tests of the 
basic control strategy and its spectroscopic verification. We ascribe to each chromophore 
a single internal harmonic vibration whose equilibrium coordinate value is displaced in 
the monomer's e-state, but whose frequency is the S£ime for both monomers and remains 
unchanged upon electronic excitation. The nuclear Hamiltonians can then be written as 
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H,=^ + ^[(«.-<0'+*']+^. (7) 
H,=^+^[«^(ft-<i)^] + ^,. (8) 

H,=^^ + ^[(q.-dy-+(q,-df] + S,. (9) 
2m 2 '- -' 

ft? is the common vibrational frequency of both intramolecular modes; and e^. are the 
donor and acceptor site energies, respectively; is the energy of the two-exciton state 
(perhaps roughly equal to £^ + fj ); and d is the displacement of the e-state nuclear 
potential curves (of importance in the calculations performed below is the dimensionless 
quEintity S = dy]mCo/2 , the coordinate displacement d divided by 2Aq^^ = ^Jl/mO) ; 

h = l throughout). qa{b) and pa{b) are the position and momentum operators for a 
vibrational mode localized on the donor (acceptor) chromophore respectively. Extending 
the Hamiltonian (1) to multiple chromophores recovers the Frenkel exciton Hamiltonian 
widely used to study energy migration in molecular aggregates.^^'^^'^"^ 

Calculations are performed using the eigenbasis of the uncoupled (7=0) 
Hamiltonian; |«)|(v„Vj)„) denotes the state in which the electronic degrees of freedom 

are in state n , while the vibrations are in an eigenstate of with and quanta in 

the donor and acceptor, respectively. The basis is truncated beyond states having more 
than a certain total number of vibrational quanta (between 10 and 24 for the parameter 
values used here). Free evolution under the EET-coupled Hamiltonian (1) is treated by 
numerical matrix diagonalization to obtain a representation of 

[t]^e-'"'=Y,\'i'j)e-"''h'j\ (10) 



in our finite basis, where j 



is an eigenket of H with energy Ej . 



We wish to predict the signals from ultrafast pump-probe difference 
measurements on the system described by Eq. (1). The measured observable is the 
contribution to the population of the one-exciton manifold (states 1 and 1') that is 
simultaneously linear in the intensities (electric field squared) of a nonresonant control 
pulse (denoted by P), an electronically resonant pump pulse (A), and a resonant probe 
pulse (Q; see Paper 1 for more details. The electric field of the pulse / has the form 
specified by Eq. (5) of Paper 1, with a Gaussian envelope function 



The first- and second-order interactions of the system with each nonzero-duration 
laser pulse are described by pulse-propagators (see Eqs. (12) and (17) of Paper 1, and 
further below). These are matrices — whose elements need be calculated only once — 
which describe a pulse's effect on the state of the system irrespective of the variable 
arrival time tj . Calculation of the pulse-propagators is simplified by assuming that the 

pulse duration Cj is short compared to the characteristic timescale 2;r/ J for excitation 

transfer, so that EET can be ignored in accounting for the action of a pulse. Using the 
pulse shape (11) and phase function (12), we find, for example. 



Mt) = exv{-(t-t,)y2c7j} 



(11) 



and a phase function 



(12) 



(K,n)jpr(-)|(v„,v, 



(13) 
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in the notation of Eq. (14) in Paper \; 5^^ is a Kronecker delta.^^ Other matrix elements 



VhVb 



of the first-order pulse propagators can be similarly obtained. The Franck-Condon 
overlap appearing on the RHS of Eq. (13) can be calculated straightforwardly using 
displacement-operator methods.^^ 

The second-order action of a pulse (such as the control pulse) includes 
contributions from matrix elements of the t5^e 

^ ' ' \ ^ ) L L (14) 

Using the eigenenergies of ground and one-exciton vibronic states, along with the 
Gaussian pulse envelope (1 1) and the phase function (12), and making a sequence of 
integration-variable changes, Eq. (14) can be evaluated analytically to give 

xexp{-cT,^A^-^K-v,f } [1 + erf (-/C7,A)] , 

where 

( y +y — \ 

^ = a^+(0 ^ "--V^ -£,. (16) 

V 2 j 

Appendix B (online)^^ details the derivation of Eq. (15) and other elements of the second- 
order pulse propagators. Possession of a compact expression for the second-order action 
of a pulse propagator — whose identification is a key practical result of this paper — 
greatly facilitates the calculation of short-pulse nonlinear optical signals by eliminating 
the need for repeated numerical evaluation of nested time-integrals like those in Eq. (14). 
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Similar expressions can also be derived for linearly chirped Gaussian pulses, whose 
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phase functions are quadratic in time. Although results of this kind might appear to be 
limited to displaced harmonic systems, they can also be applied to arbitrary anharmonic 
potentials by replacing the latter with harmonic approximations whose value, slope, and 
curvature match those of the actual potential at the instantaneous location of the wave 
packet on which the propagators act.^^ 

2. CALCULATIONS 

A. Model system with moderate electronic-vibrational coupling 

We first consider a model system with moderate vibrational displacements, 

5 = yfL5 , in the site-excited states; a value J = 0.2 co for the energy-transfer parameter; 
£ind perpendicular donor and acceptor electronic trcinsition-dipole moments. The 
monomers are assigned equal site energies, £^=£y=£ , and f 2 = + = 2f .'^^ A 

numerical search over pulse duration and center frequency, subject to a constraint of less 
thcin 5% electronic excitation probability, found that among transform-limited Gaussian 
control pulses the combination C7p = 0.140 (2;r/<») and Q.p =e-l.9l(0 generates the 

largest amplitude of nuclear motion, 0.338J = l.06SAq^^ , in the electronic ground state. 

All subsequent calculations ignore the small contribution due to control-induced 
electronic excitation. In order that the control pulse selectively excite vibrational motion 
in the "acceptor" chromophore, while the pump and probe pulses excite and probe the 
"donor," we choose a VHH polarization combination in which the control pulse is 
perpendicular to the mutually parallel pump and probe. The pump pulse, with duration 
(7^=0.1 (Itt I Q)) and center frequency Q.^=e + d^ (0,is locally resonant with the 1 <— 
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transition along the line ^„ = and with 1' <— along = . The probe pulse, with 
£7(. = 0. 25 (2;r / w) and Q.(,=e-2>5^ Q), '\s resoneint with 1 — > along the line = 2d 
and with 1' — > along the line q^=2d . These lines are the loci of vibrational outer 

turning points in the 1 and V states, respectively. 

Figure 1 shows the fractional population of the one-exciton manifold residing in 
the donor state (state- 1) of an oriented complex (one site transition-dipole vertical, the 
other horizontal) following A-pulse excitation to that state from the electronic ground 
state without and with prior impulsive-Raman excitation of coherent vibrational motion 
in the acceptor chromophore at time tp=t^- 0.56(2;r / O)) (chosen so that the ground- 
state wave packet is at its inner turning point, with = -0.338<i , when the A pulse 
arrives). Despite the fact that the acceptor-mode displacement achieved by the impulsive 
Raman process is smaller than the ideal value {qi,) = -d envisaged in Fig. 1 of Paper 1, 

the control pulse evidently slows excitation transfer significantly during the first seven 
vibrational periods. 
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FIG. 1. Relative population of donor state versus time (in vibrational periods) after electronic excitation 
from vibrational ground state (solid black curve) and from impulsively displaced vibrational wave packet 
(soUd gray curve). Corresponding dashed curves are semi-analytical predictions of weak EET-coupling 
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theory derived in Appendix A, which depend only on vibronic state populations and miss effects due to 
vibrational coherence. 

Figure 2 compares the calculated pump-probe and pump-probe difference signals 
for this system in an isotropic sample. In the absence of influence by control-pulse 
vibrational stimulation on the time-course of excitation transfer, these two signals would 
be proportional to each other at all times. The lack of simple proportionality between the 
two curves in Fig. 2 therefore demonstrates coherent control of EET. Both the survival 
probabilities and the pump-probe signals are influenced by vibrational coherences. The 
latter spectroscopic signals, though, are additionally complicated because the signal 
selects for nuclear probability density in a specified probe- window region of vibrational 
configuration space. Thus the pump-probe signals provide detailed views of the nuclear 
and electronic dynamics in this EET system, as they reflect the instantaneous spatial 
probability density in the probe-window region, rather than just the total population of a 
specific electronic state. '^^''^^ In this model system with moderately strong electronic- 
vibrational coupling, the envelope of the probe-probe signal and the pump-probe 
difference signal nevertheless roughly resemble the corresponding population traces of 
Fig. 1."*^ In particular, the envelope of the pump-probe difference signal decreases less 
precipitously during the first ~ 5(2;r/ CO) than the envelope of the pump-probe signal 
without a control pulse. 
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FIG. 2. Pump-probe signal (black curve, left ordinate in units of E^E^m jl ) and pump-probe 

difference signal (gray curve, right ordinate in units of EpE^E^fn' jl^ ) as a function of delay time 

between pump and probe pulses. Signals are calculated for an isotropic sample of the model energy-transfer 
complex with equal-energy monomers having perpendicular transition moments. 

Because the probe-pulse center frequency is strongly red-shifted from the pump 

{Q.^-Q.j^ =-4S^C0 ), both signals in Fig. 2 are dominated by their stimulated-emission 

contributions. Those contributions (not shown separately) — specified for the pump-probe 
difference by Eq. (42) and Appendix C (online) of Paper 1 — are indistinguishable on the 
scale of Fig. 2 from the total signals including also excited-state absorption and ground- 
state bleaching.'''^ 

The presence of six direction cosines (see Appendix A of Paper 1) suggests that 
the VHH-polarized pump-probe difference signal from an isotropic sample should be 
dominated by those terms which would be nonzero if the EET complex were oriented 
with one transition moment along the vertical axis and the other horizontal. In Figure 3 
are plotted those portions of the stimulated-emission signals arising from the wave-packet 
overlaps that would be generated in such ideally oriented complexes.'*^ The handful of 
overlaps contributing to the signals in Fig. 3, which is plotted on the same scale as Fig. 2, 
evidently account a major portion of the pump-probe and pump-probe difference signals 
from the isotropic sample. 
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FIG. 3. Stimulated-emission contributions to the pump-probe and pump-probe difference signals shown in 
Fig. 2 from those wave-packet overlaps that would take nonzero values if one monomer's transition 
moment were oriented vertically (parallel to the control-pulse polarization) and the other's were oriented 
horizontally (parallel to the pump and probe polarization). 

Although the main focus here is on the energy-transfer dynamics of putatively 

isolated systems, some indication of the effects of medium-induced dephasing on the 

calculated signals can be gained by introducing electronic inhomogeneous broadening.^'' 

To this end, we average the signals over a collection of isotropic dimer systems in which 

monomer site energies, £^ and e^^ , are independently distributed about a common mean 

value e . The less rapid decrease in VHH pump-probe difference signal compared to the 
HH pump-probe signal seen in Fig. 2 remains evident in the inhomogeneously broadened 
signals shown in the upper panel of Figure 4.^^ The corresponding perpendicular-probe 
cases, VHV and HV are shown there as well, and the resulting signal anisotropics — 
(VHH-VHV)/(VHH+2VHV) for pump-probe difference and (HH-HV)/(HH+2HV) for 
pump-probe — are plotted in the lower panel of Fig. 4. There is a notably lower initial 
value in the control-influenced pump-probe difference anisotropy, but it does not appear 
to decay less rapidly than the pump-probe anisotropy for this system. Appendix C 
(online)"*^ addresses the interesting dependence of the initial pump-probe difference 
Einisoptropy on the duration of pulses A and C. But the findings presented there must be 
considered approximate because the theory of Paper 1 (see its Theory section) and the 
cdculations of this paper neglect the effects of temporal A-C overlap. Such pump-probe 
overlap effects are unlikely to be entirely negligible when tc-tA = 0. 
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FIG. 4. Upper panel: Average signals from a collection of 1000 model-system dimers {6 = Vz5 , 
J = Q.2(o), whose site energies are chosen from independent normal distributions (e^ = e^- = e) with 
FWHM equal to the vibrational frequency. Pulse parameters are the same as in Fig. 2, except the temporal 
widths of the pump and probe pulses are <J^ = 0.25 (2;t / (o) and (7^ =0.5 (2;^ / (o) , respectively. Pump- 
probe signals are shown for HH and HV polarizations; pump-probe difference signals for VHH and VHV 
polarizations. Lower panel: Signal anisotropies for the inhomogeneously broadened model system, in black 
for pump-probe and gray for pump-probe difference. 

B. Dithia-anthracenophane 

Next we examine excitation-transfer dynamics and pump-probe signals from an 
EET complex with parameter values appropriate to dithia-anthracenophane (DTA), a 
covalent anthracene dimer compound that has been studied by Yamazaki and co-workers 
using time-resolved pump-probe anisotropy measurements. "^^'"^^ Two vibrational modes of 
the anthracene monomer, mode-12 at co^j^jlTtc = 385cm^' and mode-6 at 

<Wg/2;rc = 1400 cm"' , dominate the vibronic progression of its 5j <— fluorescence- 
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excitation spectrum within 1500 cm ^ of its origin at 27695 cm ' From the jet- 
spectroscopy data of Lambert et al. we abstract a dimensionless displacement 

1 for the lower-frequency vibration. We use an effective energy-transfer 
coupling constant J = (33.6cm"^ / 2) exp(^j2) = 0.0595 gj^j = 22.9 cm"^ obtained by 

dividing hcilf the population-oscillation frequency reported in Ref. 49 by the estimated 
Franck-Condon overlap between two mode- 12 wave packets separated by their relative 

displacement (-v/2J,2) in the two site-excited potential wells.^^ The two anthracene- 
monomer transition dipole moments in DTA are known to be approximately, but not 

48 

exactly, orthogonal; for simplicity, their relative angle is set to a = 7t II in our 
calculations. 

In the calculations shown here, we have taken chosen pulse widths and center 
frequencies so that only superpositions of mode- 12 vibronic states are prepared upon 
electronic excitation and only mode- 12 is susceptible to impulsive-Raman excitation in 
the electronic ground state. In particular, mode-6 is quiescent under the action of both the 
control pulse ( (7p = Q.llSilTC / (x\^ = 19.5 fs; Q.p= f -1.536^2 ) the vertically 
resonant pump ( <7^ = 0. l(2;r / 6)^2 ) - 8-66 fs ; = f + ^2 )• We therefore report 
calculations on a simplified model, referred to as DTA- 12, in which only this single 
Frcinck-Condon active mode is treated explicitly. The probe pulse is longer than the pump 
and resonant at the outer turning line of nucleeir motion in either site-excited state of 
DTA-12,with £7^ =0.5(2;^/6;|2) = 43.3fs and = f - 3^i2 ft;i2 • A VHH polarization 
scheme is adopted (as explained in subsection 2A). The control pulse generates mode- 12 
nuclear motion in the electronic ground state with amplitude 0.54<ij2 = O.SOAq^^ ■ The 
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delay between the control and pump pulses, 0.56(2;r/ a^^) ~ 48.5 fs , is chosen so that 

this control-induced wave packet reaches its inner turning point, (^j) = -0.54Jj2 , at the 

arrival time of the pump pulse. 

Calculated pump-probe and pump-probe difference signals from an isotropic 
DTA-12 sample are shown in the upper panel of Figure 5. Due to the small value, 

both excited-state absorption and ground-state bleach, along with stimulated emission, 
contribute significantly to the calculated signals. One-exciton emission signals very 
nearly equivalent to the stimulated emission components'^ of the pump-probe and pump- 
probe difference signals (Fig. 5, lower panel) could however be isolated in fluorescence 
up-conversion measurements.'^ 
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FIG. 5. Pump-probe signal (black curve, left ordinate) and pump-probe difference signal (gray curve, right 
ordinate) as function of the delay between pump and probe pulses. Upper panel shows pump-probe and 
pump-probe difference signals from simplified model of dithia-anthracenophane treating only the dynamics 
of the low-frequency mode-12. Lower panel plots stimulated-emission contributions from DTA-12. 

Despite their general similarity in form, the pump-probe difference and 
corresponding pump-probe signals in Fig. 5 are not constant multiples of each other, so 
the nuclear motion induced by the control pulse evidently affects to some degree the 
wave-packet surface-crossing dynamics in the one-exciton manifold. We turn to the 
population dynamics from an oriented DTA-12 complex to gain further insight. Figure 6 
plots the donor-state survival probability following the action of the control and pump 
pulses on a dimer with one site transition-dipole aligned with the control-pulse 
polarization (vertical) and the other aligned with the pump (horizontal). The control pulse 



induces small-amplitude motion in the acceptor — 0.54<ij2 = 0.54S^j2/mco^ 

= 0.60/ ^jlmo\2 , less than the rms width of the wave packet — so the population dynamics 
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in the presence of the control pulse differs only slightly from that in its absence. 
Moreover, the small-amplitude vibration accelerates short-time electronic excitation 
transfer in DTA-12 rather than forestalling it, by increasing the amplitude of the donor- 
to-acceptor population oscillation. As explained further in Appendix A, this behavior 
results from the fact that the vibrational displacement in mode- 12 is not large enough to 
remove the donor-state wave packet from the region of intersection between the donor- 
excited and acceptor-excited potential energy surfaces. The quasi-classical argument 
suggesting that acceptor-mode displacement should impede short-time EET (see Fig. 1 of 
Paper 1 and the accompanying explanation) does not apply in this situation. 
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FIG. 6. Donor-excited-state population dynamics for DTA-12 without (solid black curve) and with (soUd 
gray curve) prior impulsive excitation of mode-12 vibration in acceptor chromophore. Transition dipole 
moment of acceptor chromophore is aligned with vertical control-pulse polarization. Transition dipole 
moment of donor is aligned with horizontal pump-pulse polarization. Dashed curves show corresponding 
predictions of weak EET-coupling approximation (see Appendix A). 

In an effort to model the effects of electronic inhomogeneous broadening in DTA- 
12, we can introduce variable monomer site energies (as described in subsection 2A). 
Figure 7 shows the pump-probe and pump-probe difference anisotropics obtained from 
the inhomogeneously broadened stimulated-emission components. The calculated signals 
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from which these anisotropies were obtained were averaged over 1200 DTA-12 dimers 
whose site energies were chosen from independent Gaussian distributions with 100-cm' 
FWHM. Pulse parameters are the same as in Fig. 5, except <T^ = 0.25 {Itt I (O^^) = 21.6 : 

The slight control-induced acceleration of EET seen in the lower panel of Fig. 5 is 
discernible also in the shorter-period anisotropy oscillation from the pump-probe 
difference than the pump-probe signal in Fig. 7. 
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FIG. 7 Anisotropies calculated from inhomogeneously broadened stimulated emission component of the 
pump-probe signal (black) and pump-probe difference signal (gray). 

C. Downhill energy transfer 

The calculations of subsections 2A and 2B focus on homodimers, in which the 
site excitation energies of the participating monomers are the same. We next investigate 
the vibrational influence on EET in a heterodimer whose excitation donor has a higher 
site energy than the acceptor. In order to facilitate comparison with the calculations of 
subsection 2A, we use the same values of all system parameters except =8^- Icod^ . 
This difference in site energies is chosen so that the state- 1' potential energy surface 
passes through the minimum of the state- 1 potential located at iqa,<lb) = (d,0) . The 
diagonal intersection line between donor-excited and acceptor-excited potentials also 
passes through the point (q„,qi,) = {0,-d) corresponding to "ideal" acceptor-mode 
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displacement. In this "downhill" model, it is expected that electronic excitation of the 
donor chromophore at an instant when the coherently vibrating acceptor mode is at its 
inner turning point should lead to accelerated short-time energy transfer. 

The control pulse used to effect this coherent impulsive excitation has the same 
vertical polarization and pulse duration (<7p = 0. 140 (2;r/ w) ) as in subsection 2A. Its 

center frequency is down-shifted by 2(od^ to Clp = e■^, -\.9\Q) = £■^-6.9\(0 in keeping 

with the lower acceptor site energy, and the acceptor-mode wave packet is therefore 
identical to that considered previously. The horizontally polarized pump pulse has the 
same center frequency {Q.^=£^ + d^ a), duration ( <7^ = 0. 1 {l7t I (O) ), and arrival time 
relative to the control pulse as before; it is locally resonant with 1 <— along ^„ = and 
with 1' <— along q^=-d , which is largely unpopulated by nuclear probability 

amplitude belonging to either the vibrational ground state or the impulsively excited 
vibrational wave packet in the electronic ground state (given the 0.338 J range of 
control-induced vibrational motion). 

Calculations of the donor-state survival probability for a downhill complex with a 
vertically (horizontally) aligned acceptor (donor) transition moment are shown in Figure 
8. The results bear out the expectation that the initial acceptor-mode displacement should 
accelerate short-time excitation transfer and accentuate oscillatory forward-and-backward 
EET (of period ~ 6{27r I (O) ). 
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FIG. 8. Donor-state survival probabilities for oriented downhill EET model without (solid black curves) 
and with (solid gray curves) prior impulsive excitation of acceptor-mode vibration. Coherent vibrational 
excitation of acceptor is seen to accelerate short-time EET in this system. Dashed curves give the 
corresponding predictions under weak EET-coupUng approximation. 

For the downhill system, we investigate pump-probe and pump-probe difference 

signals using probe pulses of two different kinds. In order selectively to detect primarily 
the nuclear probability density in the donor-excited state, we use a horizontally polarized 
probe pulse of duration a^, = 0.25 (2;r / (O) and center frequency Q.^ =e^- 3S^ O) , which 
is resonant for 1 — > stimulated emission along the line = 2d and with 1' ^ 
stimulated emission along the line qi,=d .^^ To preferentially detect the nuclear 
probability density in the acceptor-excited state, we use instead a vertically polarized 
probe, with duration a^. = 0.25 (2;r / (O) and frequency Cl^ = - 35^ O) . This pulse is 
locally resonant with 1' — > stimulated emission near the line = 2d .^^ 

The pump-probe and pump-probe difference signals from an isotropic sample of 
the downhill dimer in the donor-detection case £ire plotted in Figure 9. The absence of 
even rough proportionality between the two signals attests to the effective influence of 
the control pulse on the time-course of short-time excitation transfer. A discernible 
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correspondence exists between the envelope of each of these signals and the 
corresponding donor-population trace in Fig. 8. But the hastening of excitation transfer 
during the first several vibrational periods by the action of the control pulse is not 
immediately obvious from the pump-probe difference signal. 
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FIG. 9. Pump-probe signal (black curve) and pump-probe difference signal (gray curve) for isotropic 
sample of downhill EET dimer using VHH polarization scheme and probe-pulse center frequency which 
select primarily for nuclear probability density in the donor-excited state. 

Interpretation of the pump-probe and pump-probe difference signals is 

complicated in the donor-detection case of Fig. 9 by the presence of significant excited- 
state absorption and some ground-state bleach, in addition to stimulated emission. These 
three contributions to both signals are shown separately in Figure 10. Stimulated 
emission (top panel) makes the largest contribution to both the pump-probe and pump- 
probe difference signals, and in both instances has a form similar to the overall signal. 
The excited-state absorption contribution to the pump-probe (pump-probe difference) 
signal has a maximum amplitude about 14% (49%) as large as stimulated emission, and 
is rather different in form. The ground-state bleach contribution to the pump-probe signal 
(bottom panel) is negligible, that to the pump-probe difference signal largely so. The 
relative sizes of these three contributions to the signals are consistent with the anticipated 
presence or absence of nuclear probability density in the vicinity of the line of local 
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resonance on the potential energy surface or surfaces from which each transition 
originates. 
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FIG. 10. Contributions from stimulated-emission (top panel), excited- state absorption (middle), and 
ground-state bleach (bottom) to pump-probe (black curves) and pump-probe difference (gray curves) to 
signals shown in Fig. 9. 

The calculated pump-probe and pump-probe difference signals from an isotropic 
sample of the downhill dimer in the alternative, acceptor-detection case are shown in 
Figure 1 1 . The strong influence of control-pulse-induced nuclear motion on the 
excitation-transfer process is again manifested by a lack of proportionality between the 
two signals. In this case, the envelopes of the two signals are very similar in form to the 
corresponding acceptor-state population (one minus the donor-state population plotted in 
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Fig. 8). Comparison of the stimulated-emission, excited- state absorption, and ground- 
state bleach contributions (not shown) shows that both the pump-probe and the pump- 
probe difference are entirely dominated by the first of these, as is consistent with 



arguments based on wave-packet access to the lines of local resonance. 
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FIG. 11. Same as Fig. 9, but with VHV polarization and probe-pulse center frequency selecting for nuclear 
probability density in the acceptor-excited state. 



3. DISCUSSION AND FUTURE PROSPECTS 

The calculations of donor-state survival probability, shown in Section 2 for three 
different EET dimers, demonstrate that impulsive excitation of coherent vibrational 
motion in the acceptor chromophore can significantly alter the time-course of population 
transfer between donor- and acceptor-excited states. The pronounced effect of an 
externally driven acceptor vibration in impeding (or accelerating) short-time EET in the 
presence of moderately strong electronic-vibrational coupling in an equal-energy (or 
downhill) dimer is demonstrated in subsection 2A (2C). The diminished effect of a 
coherent vibration on electronic population transfer under weak electronic-vibrational 
coupling is illustrated for DTA-12 by the results of subsection 2B. These illustrative 
results provide the motivation for our more extensive investigation of the various 
corresponding spectroscopic pump-probe difference signals. 
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The at-best indirect correspondence between the nominally observable total 
population of the donor-excited state and the spectroscopic signals that would actually be 
measured in an ultrafast pump-probe or pump-probe difference experiment suggest that 
detailed system- specific calculations will generally be required to assist in the 
interpretation of future experiments along the lines proposed here and in Paper 1. 

On the other hand, the calculations of subsection 2A demonstrate that the pump- 
probe difference signal (for appropriately chosen pulse center frequencies) can be 
dominated by its stimulated-emission component and that this signal component from an 
isotropic sample does not differ greatly from that for an oriented complex (which 
comprises just two wave-packet overlaps"^^). These findings suggest that similar 
simplifications should be possible in the calculation cind analysis of the full nl-WPI 
difference signals for which the general expressions are given in the preceding paper. 
Domination of nl-WPI difference signals from EET systems by a small number of 
overlaps should greatly facilitate their interpretation in terms of the amplitude-level 
dynamics of surface-crossing wave packets. 

The signal calculations presented here focus mainly on isolated model systems 
most closely representing energy-transfer dimers in the gas phase. Neglect of rotational 
motion is reasonable given the nanosecond-to-picosecond rotational periods estimated for 
molecules such as dithia-anthracenophane.^'^'^^ But following the strategy of Ref. 10, the 
effects of medium-induced inhomogeneous broadening were mimicked in some cases by 
averaging calculations from EET systems having variable independent monomer site 
energies selected from a Gaussian distribution. Dynamical effects of embedding EET 
dimers in a liquid or crystalline environment could be simulated by combining Redfield^^' 
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' ' ' ' or other relaxation theories with expUcit dynamical treatment of strongly 
driven nuclear degrees of freedom, or perhaps by semiclassical wave-packet methods.^"^ 
The full nl-WPI signals considered in Paper 1 depend on the preservation of electronic 
coherence over the relevant intrapulse-pair delays. But pump-probe difference 
measurements of impulsive vibrational control should be primarily sensitive to 
vibrational dephasing, and should therefore be applicable to condensed-phase systems on 
a timescale of tens of picoseconds. 

The control strategy investigated here is motivated by the semiclassical notion 
that a Gaussian wave packet displaced from the - line of intersection between 

donor- and acceptor-excited electronic states of equal site energy should be less 
susceptible to energy transfer than one sitting astride that line. As seen in subsections 2A, 
this strategy is effective under fairly strong electronic-vibrational coupling. On the other 
hand, subsection 2B shows that this simple displacement method is ineffective under 
weak electronic-vibrational coupling, where the spatial width of the wave packet is 
smaller than the excited- state shift in the equilibrium value of the vibrational coordinate, 
and the semiclassical picture of EET breaks down. 

This failure under weak electronic-vibrational coupling does not, however, 
preclude the exertion of vibrational control over excitation transfer in such systems. With 
reference to Appendix A, we note, for instance, that the Franck-Condon overlap for 
DTA-12, ((2^)j, I (2^X) = -0.03506 , happens to take a very small value. This small 

overlap should lead to slow excitation transfer from any state of the form 1 1) | with 

I ) = f/, ® /, I (0, 2 J, ) = ® I (2„ J, > - ^ I (1, 1, ) + i I (0„ 2,\)} The donor- 

state survival probability following the preparation of such an initial state in DTA-12 is 
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plotted in Figure 12. For the example shown, we have imagined displacing the wave 
packet I i^o) = I (Os 2a)o) by -d along prior to short-pulse absorption to the donor- 
excited state, so that \yfi) = e'""'^ \ (O^l^)^) = e'^^''''^ \ (0^ 2^)i) . Although this illustration 
defers the question of how the initial state |0)|(0^ '^a)q) could be prepared in practice, the 

greatly diminished initial decay of the donor-state population observed in Fig. 12 
suggests that effective vibrational control over EET might well be possible even in the 

presence of weak electronic-vibrational coupling. 

qjd 

-4 -2 2 4 
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FIG. 12. Illustration of possible vibrational-control strategy for the DTA-12 model with its weak electronic- 
vibrational coupling. Upper panel shows spatially translated = 2 vibrational wave function at instant of 
transfer to donor-excited state. Lower panel compares subsequent donor-state survival probability (gray 
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dashed line) to that following direct excitation from vibrational ground state (solid black line). 
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APPENDIX A: POPULATION TRANSFER FOR J«(o 

Oscillations at the vibrational frequency are generally less prominent in the 
survival probability traces of Figs. 1, 6, and 8 than they are in the corresponding pump- 
probe and pump-probe difference signals of Figs. 2, 5, and 9, respectively. While the 
latter signals are sensitive to the nuclear probability density in localized spatial windows, 
the former depend on the nuclear dynamics only through its effect on the detailed time- 
course of population transfer between the site-excited electronic states. When the EET 
coupling constant J is small compared to the vibrational energy spacing, energy transfer 
occurs predominantly between vibronic states of nearly equal energy in the two site-state 
meinifolds, if such energy matches exist. In the absence of interfering treinsition 
amplitudes from multiple vibronic levels of the donor, for instance, vibrational 
coherences do not impress themselves on the population evolution of a given vibronic 
level in the acceptor. The total acceptor population and the survival probability of the 
donor excitation are then largely devoid of vibrational-frequency oscillations. 

Let us consider the electronic population dynamics in the one-exciton manifold of 
a system governed by 

^.„. = |l)^i(l| + |l>r(l1 + ^ , (17) 

in which H\ and Hr are given by Eqs. (7) and (8), respectively, with < , and 
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(18) 



The population evolution can be efficiently described by switching to symmetric and 
antisymmetric combinations of the donor and acceptor nuclear coordinates and momenta: 



(see Figure 13). In terms of these operators, 



'1 



2m 



+ 



and 



P + P 



2m 



+ - 



mco 



+ 



Q 

' V2 



(19) 



(20) 



We denote by I (v^V^)i) (or I (V^V^)i') ) an eigenket of (or Z/^- ) with energy 



aXy^ +v^+l) + e^ (or aXy^ + v^+V) + e^,). Note that the equilibrium values of (Q^ , ) 

are (d I -spl, -d I yfl) and (d I yfl, d I V2) in state- 1 and state- 1', respectively; only the 

equilibrium value of the antisymmetric coordinate differs in the two states. 
We wish to determine the survival probability 

p,(0 = (?^,|(i|e^^-'|i)(i|e-'"-'|i>|?^i> 



for an initial state in the donor-excited manifold (with i\i/^\Yi) = '^ )■ This can be 



accomplished by solving the equation of motion 



dt 



for the one-exciton time-evolution operator in the interaction picture, where 

V(t) = e'^»one-y)ty^-i(H,„-y)t _ 



(21) 



(22) 



(23) 
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The interaction-picture "Hamiltonian" (23) has matrix elements 

{{yyA)A{A^it)\l)\{v,vA) = Je'''-^^^^^^^^^^ (24) 

and 

((v,vjJ(i|y(0|i)|(vXX') = -/^''^^'-^''^'""'^^-^^''((^^ • (25) 

Elements (24) and (25) both oscillate at frequency (fj. -£]) + 0)(y\ - v^) ; in the weak- 
coupling case (7 « O)) , they are therefore ineffectual unless e^. + v'^^O) ~ v^co . 

In order to investigate the simplification that becomes possible when the site 
energies are equal or differ by multiple vibrational quanta, we assume -e^^ = Net), hi 

this situation, it is reasonable to omit the time-dependent elements of V , and replace it by 

= 7 X 1 1 ' + A^)r > ((^a + A^)i' | (^a )i > ((v. > ^a )i | (l | + H.c. , (26) 

which is time-independent. Substitution in Eq. (22), integration, and return to the 
Schrodinger picture then yield 

Under this approximation, the survival probability (21) becomes 
Pi (0 = {y/, I (1 1 e'" 1 1) e'"^'e-'"^' (l | e"'*' \\)\y/,) 

= ({^^1 1 ^1 1 cos i)? 1 1) (l I cos i)? 1 1) I j^^i ) 

= ({^^i|(l|cos^i)?|l)|{^^i) 

^■ + ^ E \{{ysyA)Mt^os[ltJ{{v^ + N\\iv,\)'] . (28) 



2 2 



VS'^A 



We want to evaluate Eq. (28) in situations where, as a result of impulsive control- 
pulse excitation, the initial wave packet takes the approximate form of the lowest state-0 
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vibrational eigenket displaced a distance / along the acceptor mode (as illustrated in 



Fig. 13): 



|?^,) = e-""'^|(OA)o) 



(29) 



We also have 



|(v.vJ,) = e-"''|(v,vJo) 



(30) 



so the two-dimensional overlap in Eq. (28) becomes 



4. 



mO) 



if-d) 



mO) 



if + d) 



(31) 



The survival probability is thereby obtained in the form 



1 1 -'^(fW^ 1 



2 2 



mO) 



if + df 



cos[2tJ{(v^+N),\(vJ,)] . (32) 



The requisite one-dimensional Franck-Condon overlaps are straightforwardly evaluated 
using the relation 



((V, + N), I (vj,) = (-V2 Sfe-'' y i-2S'Y ^^""^ ^ -l^^ . (33) 
\v A n\yAH/ ' A;!(A; + iV)!(v^-it)! 
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FIG. 13. Symmetric and antisymmetric vibrational coordinates plotted on donor and acceptor coordinate 
axes. The lowest-energy vibrational state in the donor-excited [acceptor-excited] electronic state is centered 

at {q^ ,q^^) = {d,Q) [ (^^ , ) = (0, rf) ]. Also shown are the location (at the origin) of the lowest 

vibrational state in the electronic ground state and the location of a ground- state wave packet displaced by 

/ along the -axis . 

In the case of "ideal" displacement f =-d , the approximate survival probability 
(32) reduces to 

Pi(0 = ^ + |cos[2fy((iVJ,|(0jJ] , (34) 

and population transfer proceeds as if the state- 1 wave packet were the lowest vibrational 
level of the donor-excited electronic state. This behavior is to be compared with the 
survival probability following direct Franck-Condon excitation (/ = 0) , 

P,(0 = ^ + ^^'-'^/^|;^^^cos[2f7((v,-hA^),|(vJ,)] , (35) 
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in which non-zero populations in higher donor-excited levels also contribute. The relative 
speed of short-time EET evidently hinges on the size of the overlaps (33) for contributing 
donor-excited vibronic states. 

The weak electronic-vibrational coupling limit ( « 1 ) is exemplified for the 
case of equal site energies by the DTA-12 calculations of subsection 2B. In this system, 
the vibrationally excited antisymmetric-mode overlaps are all smaller in size than 
((Oyi)i' I (Oyi)i) ' the short-time donor-population falls off more rapidly following 
acceptor-mode displacement (Eq. (32) with / = -0.54 J ) than it does under Franck- 
Condon excitation (Eq. (35)). These predictions of the weak EET-coupling 
approximation are consistent with the rigorously calculated population dynamics shown 
in Fig. 6. In addition, they differ only slightly from the survival probability calculated by 
means of Eq. (28) using the actually generated by the sequential action of the 

control and pump pulses, rather than the approximate form (29) (results not shown). 

The model system studied in subsection 2A features equal site energies, moderate 
electronic-vibrational coupling (S^ = 2.5) , and moderate energy-transfer coupling 
(J = 0.2(0) . For this system under "ideal" impulsive displacement of the acceptor mode 
(/ = -d) , Eq. (34) with A^^ = predicts 



1 1 f2jA 

2 2 \e J 



(36) 



According to Eq. (35), vertical donor excitation without prior vibrational displacement 
leads to 
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(37) 



Here = 1 - 4 all have larger overlaps than = , so the donor-state population falls 
off more rapidly after direct Franck-Condon excitation than it does following prior 
acceptor-mode displacement. For the displacement, / = -0.338<i , achieved by the 
impulsive Raman process described in subsection 2A, Eq. (32) predicts 
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+ 0.5478COS 
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(38) 



This function is plotted along with the rigorously calculated population in Fig. 1; a 
similar form that agrees slightly better with the rigorous result (not shown) is obtained 
from Eq. (28) evduated with the numerically calculated components ((VjV^)i | Y\) 

generated by the actual impulsive Raman process and subsequent pump-pulse excitation. 

In the intermediate electronic- vibrational coupling case ^ = 1 , the Franck-Condon 
overlaps for A'^ = given by Eq. (33) depend only weakly (in magnitude) on the 
antisymmetric-mode quantum number, taking values \l e , -\l e, -II e, -l/3e , and 
1 / 3e for v^=0 through 4, respectively. These overlaps give rise to survival probabilities 



1 1 (2Jt^ 
R(t) = - + -cos 
'22 



(39) 
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for / = , so the acceptor-mode displacement has only a small effect in this case. This 
pattern suggests that simple displacement of a moderately Franck-Condon active acceptor 
vibration, such as mode-6 of anthracene (with = 1.05 <»g = 1400cm~' in DTA, 

would be largely ineffective as a means of controlling electronic excitation transfer. We 
have verified this prediction with numerical calculations (not shown) of the donor-state 
survival probability and pump-probe signal for a DTA-6 model analogous to the DTA- 12 
results of subsection 2B. 

In the downhill model considered in subsection 2C, the relevant overlaps are 
{(Va +5)r Kv^X) = -0.4189, -0.1710, 0.1371, 0.2706, and 0.2037 for Va = 0-4, 

respectively. The corresponding population-weights to be used in Eq. (32) with / = 
and / = -0.338<i are the same as those appearing in Eqs. (37) and (38), respectively. 
The resulting, approximate predictions for (t) axe plotted in Fig. 8 along with the 

rigorously calculated traces. 
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APPENDIX B: VIBRONIC MATRIX ELEMENTS FOR SECOND-ORDER 
PULSE PROPAGATORS 

In going to the pump-probe limit of a nl-WPI experiment, the interpulse delays for pulses 
within a phase-related pair arc set to zero. The treatment described in Paper 1 (see Eq. 
(13)) properly accounts for the non-negligible effects of pulse overlap that enter in this limit. 
We construct second-order pulse propagators for ground-state bleach (GSB), stimulated 
emission (SE), or excited-state absorption (ESA) depending on the particular action of the 
pulse (0— >£— >0,£— >0— >£,or£— >2— >£ respectively). 

Propagators for GSB take the form 

pT\oo- t,)pr\t,- to = j dt, I cit,/,(t,)//(ti)e^*^(*^)-^*^(*^) 



- oo — oo 



(B.l) 



^^-iHoiti-t2)^-iH^{t2-ti)^-iHi(ti-ti)^-iHo{ti-ti) 



To calculate the matrix elements for this operator in the eigenbasis of Hq, we introduce 
completeness relations in the eigenbasis of Hi, 



(^a,'^b)o pf^\oo; t2)pf°\t2; ti) {Ua,'^b)i 

^ ^_Eim^ j j ^^^g-(ti-t,)72a2g-(t2-t,)Y2a|gin,(t2-t,)-in,(ti-tj) 

— oo — oo 

= K,J>b i^—Y^ ( (( ^a)e\ ( ^'0)9) 



00 T2 



X 



j dT2 J dri exp -Fi (a;(i/a - i/a) - Qj) Tij 



-00 —00 



X exp \ --^ - i (a;(i/a - i^a) + £1 - ^i) T2 



2aj 



(B.2) 
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A closed-form expression can be derived for nested time integrals of the form appearing in 
Eq. (B.2), due to the Gaussian nature of the integrand. The first step is to change the order 
of integration, 



OO T2 



I{a,(3)= J J dri exp i^--^ + iaTi - - i Pt2^ 



-OO —CO 
OO OO 



.2 2 
^2 



(B.3) 



J dri J dT2 exp <| + iari - ^ -if3r2 



— OO ri 



We now make a change of variables, letting t = T2 — ti 



OO OO 



f f ( + + 2 T Ti 1 

/(a,/3)= J dn J dr exp^-^^+ian- ^ i p {t + n) j , (B.4) 



-00 

and then change integration order again to obtain 



00 ^ 00 ^ 

'(q;,/3) = J drexp i^--^ - i Pt^ J dri exp + i (a - /3)ti - . (B.5) 





The inner integral is solved in the usual way, giving 





00 



I T 



(B.6) 



which can be written in terms of the error function, 

/(a, /3) = TT a] exp I - ^ (a^ + I (^1 - erf I ^ (a + /3) I ^ . (B. 7) 



We can now return to the vibronic matrix element for the GSB propagator, which by Eq. 
(B.7) can be written as 

= Kfi, (^^) ^("""^sl ( ^( ( "^"^s) (B.8) 

X / - i/a) + £1 - 0/, L0{pa - Va) + ^1 - ^l) ■ 
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Matrix elements of the operator pf^ \oo; t2)pf ti) follow analogously. A GSB propa- 
gator also describes the action of the pre-resonant control pulse, whose effect on the system 
is primarily the creation of a ground-state wave packet, and Eq. (15) of the main text is 
equivalent to Eq. (B.8). 

The ESA propagators are obtained in a similar fashion, 

2 OO t2 
— OO — OO 



(iE Tn\^ 



X / {uj{pb - Pb) + (£2 - £l) - ^I: l^ih - l^b) + (^2 - £l) - ^l) , 

(B.9) 

and the corresponding elements of pj (oo; t2)pj {t^-, ti) follow by direct analogy. Unlike 
in the case of GSB, we now have the possibility that the final and initial electronic states 
will be different. We then need matrix elements of the form 



J dt2 J (iiie~^*'"*^^V2^ie"(*2~*^)V2<^?e^"^^*2~*^^~^"^^*i~*^^ 



— OO — OO 



2 



I (i^a)e) ((i^6)e| ( i^b)^) 

X / {uj{Pb - I'b) + (^2 - £l) - <^{i^a - l^a) + (^2 " ^l') " ^l) ■ 

(B.IO) 
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Finally, the SE propagators are written 

^ ^Eim ^ j j ^^^^-(ti-tj)y2af^-it2-tir/2ai^-in,(t2-ti)+ini{ti-ti) 

— oo — oo 

(B.ll) 

and 

((i^a,i^6)l'|pi^°^(oo; t2)pf^\t2; h) |(Pa,i^&)i) 

^ l^ ^^/m ^ j j ^^^^-(t,-ti)y2al^-{t2-tif/2af^-ini(t2-tj)+in,(h-tj) 

— oo — oo 

= (^^^) I 

(B.12) 

APPENDIX C: INITIAL ANISOTROPY VALUES FOR PUMP-PROBE 
AND PUMP-PROBE DIFFERENCE SIGNALS 

For the pump-probe signal, the anisotropy is defined by 

, , HH-HV 

= HH + 2HV • (^-'^ 

Using well-known orientational factors^, we can write the stimulated emission component to 
the numerator above (assuming the monomers have orthogonal transition dipole moments) 
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as 



HH - HV =^ Re (4 ({a(10)},|{c(10)c(01)a(10)},) + 4 ({a(l'0)},,|{c(l'0)c(01')a(l'0)},,) 

- 2 ({a(l'0)},|{c(10)c(01)a(l'0)},) + 3 ({a(l'0)},,|{c(l'0)c(01)a(10)},,) 

+ 3 ({a(l'0)}J{c(10)c(OlXlO)}i) + 3 ({a(10)}i,|{c(l'0)c(01)a(l'0)}i,) 
+ 3 ({a(10)}J{c(10)c(01')a(l'0)},) - 2 ({a(10)},,|{c(l'0)c(01')a(10)},,)) 
= ^ Re (4 ({a(10)} J{c(10)c(01)a(10)},) - 2 ({a(l'O)} J{c(10)c(01)a(l'0)}i) 
+ 3 ({a(l'0)},,|{c(l'0)c(01)a(10)},,) + 3 ({a(l'0)},|{c(10)c(0l')a(10)},)) . 



(C.2) 

In the second equality, use has been made of the fact that for equal-energy homodimers the 
overlaps retain their value when the labels 1 and 1' are interchanged. The denominator is 

HH + 2HV=jRe(({a(10)}J{c(10)c(01)a(10)}i) + ({a(l'0)}J{c(10)c(01)a(l'0)}i) 
+ ({a(10)},,|{c(l'0)c(0l')a(10)}„) + ({a(l'0)}„|{c(l'0)c(Ol')a(l'0)}„)) 
=^Re(({a(10)}J{c(10)c(01)a(10)},) + ({a(l'0)},|{c(10)c(01)a(l'0)},)) . 

(C.3) 

When the interpulse delay tcA is zero, overlaps like ({a(l'0)}^|{c(10)c(01)a(l'0)}]^) vanish, 
as there is no interval of free evolution during which energy transfer can take place. This 
simplification gives a compact expression for the initial anisotropy 

rpp(0)-0.4 + 0.3 R,„„(io)j^,{,(io)e(01)a(10)},> ' '^'^^ 

which reduces to 0.4 or 0.7 when the ratio of wave-packet overlaps in the second term is or 
1, respectively. This ratio is essentially zero under the conditions of our simulations (which 
use a red-shifted probe pulse) due to the fact that the spatial region in which the probe 
pulse is resonant with the 1 — > transition (near qa — 2d) overlaps the region in which it 
is resonant with the 1' transition (g?, = 2d ) far from the location of the wave packet 
created by the pump pulse. Were the pump and probe to have the same center frequency, 
or were the probe pulse sufficiently short that its center frequency became irrelevant, the 
value of rpp(O) = 0.7 would be obtained (see Figure panel A). 
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The anisotropy for the pump-probe difference signal, in which the system first interacts 
with a vertically polarized control pulse, is given by 

, , VHH - VHV ,^ , 

= VHH + 2VHV • ^^-^^ 

Omitting those overlaps that depend on energy transfer between pulses, and again making 
use of the symmetry with respect to interchange of 1 and 1', we find 

VHH - VHV = ^ Re (2 ({a(10)}J{c(10)c(01)a(10)p(01')p(l'0)},) 

+ ({a(l'0)},,|{c(l'0)c(01)a(10MOlMlO)}i,) + ({a(l'0)}i,|{c(l'0)c(01)a(10)p(0l')p(l'0)},,) 
+2 ({a(10)p(01')p(l'0)}J{c(10)c(01)a(10)}i) + ({a(l'0)p(01)p(10)}i,|{c(l'0)c(01)a(10)}i,) 
+ ({a(l'0)p(01')p(l'0)},,|{c(l'0)c(01)a(10)},,)) 

(C.6) 

and 

VHH + 2VHV = ARe(9 ({a(10)}i|{c(10)c(01)a(10)p(01)p(10)}i) 
iUo 

- ({a(l'0)},,|{c(l'0)c(01)a(10)p(01)p(10)},,) - ({a(l'0)},,|{c(l'0)c(01)a(10)p(0l')p(l'0)},) 
+ 13 ({a(10)},|{c(10)c(01)a(10M01'Ml'0)}i) + 9 ({a(lOMOlMlO)} J{c(10)c(01)a(10)}i) 

- ({a(l'0)p(01)p(10)},,|{c(l'0)c(01)a(10)},,) - ({a(l'0)p(0l')p(l'0)},,|{c(l'0)c(01)a(10)},) 
+13 ({a(10)p(0l')p(l'0)},|{c(10)c(01)a(10)},)) . 

(C.7) 

In the case that all of the overlaps listed in Eqs. (C.6) and (C.7) are equal in value (which 
happens when the pump and probe are arbitrarily short), the initial anisotropy is 0.7. Nu- 
merically, we find an initial anisotropy of 0.69 when o"^ and uc are set equal to one hundredth 
of the vibrational period for the equal-energy dimer used in Section 2A of the main text (see 
Figure panels B and C). 

When the overlaps containing c(l'0)c(01) do not contribute, for the same reason as de- 
scribed above, we have 

rppD(O) = 

I (C.8) 

13 9 ^'^((H10)}i|M10M01)q(10)p(01)p(10)}i) + ({«(10)p(01)p(10)}i|{c(10)c(01)«(10)}i)) 
Re(({a(10)}i|{c(10)c(01)a(10)p(01')p(l'0)}i) + ({a(10K01')p(l'0)}i|{c(10)c(01)a(10)}i)) 
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In general, this expression gives an initial anisotropy that changes with pulse and molec- 
ular parameters, in contrast with the initial value of 0.4 often found in pump-probe ex- 
periments. Were the control pulse arbitrarily short, it would be ineffective in driving nu- 
clear motion in the electronic ground state, and |{a(10);7(01)p(10)}^) would be identical 
to|{a(10)p(01')p(l'0)}i). In this limit, Eq. (C.8) reduces to a value of 7/22 w 0.32 for non- 
impulsive pump and probe pulses (see Figure panel B). That this value for an impulsive 
control pulse is different from 0.4 is not unexpected. The control pulse creates a copy of 
the ground-state nuclear wave packet only in those molecules in which one of the monomers 
has a nonzero transition dipolc component along the V axis. It is this subset of the total 
isotropic population upon which a pump-probe experiment is conducted. In the situation 
of interest - where the control pulse does generate motion on the electronic ground state - 
we find no single limiting value for Eq. (C.8) (see Figure panel C). 

It should be recalled that our numerical simulations explicitly ignore overlap between 
the pump and probe pulses, and our results are therefore only strictly valid for interpulse 
delays longer than the pulse lengths. Numerical calculations that rigorously include overlap 
between the pump and probe (not included here) show that neglecting this overlap tends 
to result in an overestimation of the initial anisotropy for the pump-probe difference. For 
example, using the same pulse and molecular parameters as in Fig. 4 of the main body of 
the paper and neglecting pulse overlap gives an initial anisotropy of 0.25 (which is largely 
unaffected by inhomogeneous broadening), whereas properly treating the effects of pulse 
overlap reduces this value to 0.14. Thus the "initial anisotropics" plotted as a function of 
pump and probe pulse duration do not strictly depict the true initial anisotropy. Nonetheless, 
due to the presence of a single vibrational period in our two-mode model system, they do 
hold significance as the anisotropy when the interpulse delay is equal to the vibrational 
period, if energy transfer is neglected on this time scale. 



See, for example, A. Tokmakoff, "Orientational correlation functions and polarization selectivity for non- 
linear spectroscopy of isotropic media. I. Third order," J. Chem. Phys. 105 1 (1996); A. Tokmakoff, 
"Orientational correlation functions and polarization selectivity for nonlinear spectroscopy of isotropic me- 
dia. II. Fifth order," J. Chem. Phys. 105 13 (1996). 
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FIGURE. Initial anisotropy neglecting pulse overlap 
(or the anisotropy one vibrational period after the 
pump pulse arrives, neglecting energy transfer) as a 
function of pump and probe pulse duration for the 
equal-energy model system with moderate electron- 
vibrational coupling (5^ = 2.5). The pump pulse 
is vertically resonant at the ground-state equilibrium 
geometry, while the probe is resonant at the outer 
turning point for nuclear motion in the excited elec- 
tronic state {^A = El + ^^f^ and = Ci — S6^uj). A) 
The pump-probe anisotropy takes a value of 0.7 for 
an impulsive probe, and 0.4 otherwise. B) The pump- 
probe difference anisotropy using an impulsive control 
pulse, showing the same behavior as (A) but with a 
lower limit of 7/22 rather than 0.4. C) The pump- 
probe difference anisotropy using a control pulse opti- 
mized to generate coherent ground-state nuclear mo- 
tion (as described in the main text) shows more com- 
plicated behavior than (A) or (B). The anisotropy is 
high for an impulsive probe, however as the probe 
pulse becomes longer the initial anisotropy does not 
decrease monotonically. 
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